Ghost imaging lidar via sparsity constraints 
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For remote sensing, high-resolution imaging techniques are helpful to catch more characteristic information of 
the target. We extend pseudo-thermal light ghost imaging to the area of remote imaging and propose a ghost 
imaging lidar system. For the first time, we demonstrate experimentally that the real-space image of a target 
at about 1.0 km range with 20 mm resolution is achieved by ghost imaging via sparsity constraints (CISC) 
technique. The characters of CISC technique compared to the existing lidar systems are also discussed. © 
2013 Optical Society of America 

OCIS codes: (110.2990) Image formation theory; (280.3640) Lidar; (100.3010) Image reconstruction techniques. 
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As an important detection tool, lidar has been widely 
used in remote sensing in recent decades. Traditional 
imaging lidar, according to the measurement modes, can 
be classified into two types: scanning imaging lidar and 
non-scanning imaging lidar (e.g. flash radiography) [1-3]. 
Scanning imaging lidar obtains the real-space image of 
the target by scanning the target region point-by-point 
with a pulsed laser [2]. Therefore, it is difficult to im- 
age moving target with high-speed. Non-scanning imag- 
ing lidar [3], which is characterized by an imaging sys- 
tem with high resolution and a pulsed flash laser, can 
cover the whole field of the target and obtain the tar- 
get's real-space image in a single exposure. However, be- 
cause the light intensity reflected by the target is divided 
into many small pixels of charge-coupled device (CCD) 
camera, the detection sensitivity is low and the detec- 
tion distance of non-scanning imaging lidar is limited by 
the signal-to-noise ratio (SNR) of the received photons 
distributed on the whole imaging plane and the aperture 
of the imaging system. In addition, for scanning imaging 
lidar the imaging resolution is limited by the Rayleigh 
criterion of the emitting aperture [4] while by the numer- 
ical aperture of the receiving system for non-scanning 
imaging lidar. 

Recently, there are many researches on the remote 
sensing with ghost imaging technology [5-11]. Combined 
with sparsity constraints, super-resolution [12-14], com- 
pressive sensing [14-19], compressive radar [20] and other 
techniques (http://ecos.maths.ed.ac.uk/SPARSll/) are 
possible. In this Letter, a practical ghost imaging lidar 
via sparsity constraints (CISC lidar) system was pro- 
posed and high-resolution imaging was experimentally 
demonstrated at a distance about 1.0 km range. 

Fig. 1 presents the setup of the proposed system. The 
source, which consists of a 532 nm coherent solid-state 
pulsed laser with 10 ns pulsed width, a rotating diffuser 
(ground glass disk) and a set of field lens, forms speckle 
field at stop 2 (field stop). The light emitting from the 
source is divided by a beam splitter (BS) into an object 
and a reference paths. In the object path, the light prop- 
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Fig. 1. Experimental setup of CISC lidar. 



agates through the objective lens /o and then to the tar- 
get. The photons reflected by the target are received by a 
light concentrator (a Cassegrain telescope with the aper- 
ture 420 mm and the focal length 5 m) and then passes 
through an interference filter with 1 nm half-width into 
a photomultiplier tube (PMT). In the reference path, 
the light goes through the reference lens /i and then to 
a CCD camera. In this system. Stop 2 is placed on the 
conjugate plane of both the target and the CCD camera, 
which can control the field of view (FOV) on the target 
plane. The field lens is used to increase the utilization 
rate of light energy and generate a visual aperture stop, 
which is signed as stop 3 in Fig. 1. Therefore, the trans- 
verse size of light beam at the objective and reference 
lens is controlled by the stop 1, which also ensures that 
the entrance pupil is exactly the same for the lens /o 
and /i . In addition, the PMT is used to transform light 
signals reflected from the target into electric signals and 
the interference filter is used for inhibiting background 



light. 

By exploiting the image's sparsity constraints, CS re- 
construction techniques usually yield, as predicted by 
theoretical analysis and confirmed by experiments, bet- 
ter results when the target is sparse in the representation 
basis [14-19]. The reconstruction of GISC lidar will be 
formulated in the CS framework. By reshaping each of 
the speckle intensity distribution (to x n pixels) recorded 
by the CCD camera (Fig. 1) into a row vector (1 x N, 
N = mxn), the measurement matrix A (M x N) is ob- 
tained after M observations. Meanwhile, the intensities 
recorded by the PMT in the object path are arranged as 
a column vector Y (M x 1). If we reshape the unknown 
target {m x n pixels) into a column vector X {N x 1) and 
X can be represented as X = ^E*- a such that a is much 
sparser (^P denotes the transform operator to the sparse 
representation basis), then the target's image can be re- 
constructed by solving the following convex optimization 
program [13,18,21]: 



X = ^ • a; which minimizes : — 

2 



|Y-AX||^ 



(1) 



where r is a nonnegative parameter, ||T^||2 and ||V^||j^ de- 
note the Euclidean norm and the i!i-norm of V, respec- 
tively. 

To experimentally demonstrate the characteristics of 
GISC lidar, the concrete parameters of GISC lidar in 
the experiments are as follows: the focal length of the 
lens /o is 5 m and the magnification of the lens /i is 
1 x . The FOV of the receiving system and the emitting 
system at Iq =900 m range is about 2 m and 1 m, respec- 
tively. The emitting aperture (stop 3) in the test path 
is 18 mm (measured value), then the theoretical resolu- 
tion (Rayleigh criterion) for traditional imaging system 
mm~ 32.5 mm. According to the parameters 
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of emitting system, the pixel size of the GCD camera 
in the reference path is set as 27.6 /zm x 27.6 fim. The 
image is reconstructed using the gradient projection for 
sparse reconstruction algorithm [21]. 

The first demonstration of GISC lidar was performed 
using our designed targets mounted on a building located 
about /o'=900 m away. The targets arc highly reflective: 
a standard Chinese vehicle license plate and a set of res- 
olution panels. The line width of the characters on the 
vehicle license plate, as shown in Fig. 2(c), is about 10 
mm. The imaging area on the CCD camera is 140x72 
pixels, which means that the imaging area on the target 
plane is 695 mm x 357 mm. Fig. 2(d) presents the im- 
age of the vehicle license plate reconstructed by GISC 
lidar. The resolution panels, as shown in Fig. 2(c), are 
divided into three groups (six-slit, three-slit and double- 
slit) and their centcr-to-center separation between slits 
are 20 mm, 40 mm and 60 mm, respectively. The imaging 
area on the CCD camera is 128x92 pixels, which means 
that the imaging area on the target plane is 635 mm x 
457 mm. The reconstructed results of the resolution pan- 
els are illustrated in Fig. 2(f) and Fig. 2(g). From Fig. 
2(f) and Fig. 2(g), the resolution panels with the reso- 
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Fig. 2. Experimental reconstruction results for high- 
reflection targets we proposed at 900 m range (with 3000 
measurements) . (a, b) The original target plates imaged 
by a camera and a telescope, respectively; (c) the con- 
crete sizes of a standard Chinese vehicle license plate; (e) 
the concrete sizes of a set of resolution panels; (d) and 
(f) are the targets' images reconstructed by GISC lidar 
and the targets arc all represented in the space basis; (g) 
the cross-section of the rectangular selection box in (f ) . 



lution 20 mm, can be clearly differentiated, which also 
demonstrates that super-resolution imaging is achieved 
by our GISC hdar. 

Another demonstration of GISC lidar was conducted 
to image the overhead target of a 200 m tower located 
about /o'=720 m away. The target's images, taken by a 
camera and a telescope, are shown in Fig. 3(a) and Fig. 
3(b). Fig. 3(c) presents the reconstruction result of the 
target (where the target is represented in 2D-DCT ba- 
sis). In Fig. 3(c), due to the ultra- low reflectivity and the 
specular reflection, the three black antennas are invisi- 
ble. 

Compared with traditional imaging lidar, GISC lidar 
has both the advantages of scanning imaging and non- 
scanned imaging lidar. Similar to scanning imaging lidar, 
GISC lidar has high detection efficiency and long detec- 
tion distance since all photons collected by the concen- 
trator illuminate the same PMT. Similar to non-scanning 
imaging lidar, the laser pulse emitted form the GISC li- 
dar covers the whole detection field, therefore, the im- 




Fig. 3. Reconstruction results of a natural target (the 
overhead target of a tower) at 720 m range (with 3000 
measurements), (a) and (b) are the image of the target 
imaged by a camera and a telescope, respectively; (c) is 
the target's image reconstructed by GISC lidar and the 
target is represented in two-dimensional Discrete Cosine 
Transform (2D-DCT) basis. 



age can be reconstructed without scanning the target 
and the imaging of targets with high-speed is possible 
even when the sampling number is far fewer than the 
pixel number the desired resolution needed [17-19]. In 
addition, the CCD camera in the reference path, which 
limits the sampling speed at present, can be omitted 
by using the techniques such as computational ghost 
imaging [22, 23] and encoding the pseudo-thermal light 
source [24] , and the approaches such as image separation 
reconstructions [25,26] are also helpful to the imaging of 
targets with high-speed. Furthermore, as indicated by 
the results in Fig. 2 and in Ref. [13], GISC lidar can also 
realize supper-resolution imaging. 

In conclusion, we experimentally demonstrate a novel 
imaging lidar by combining GI method with sparse and 
redundant representations. We show that GISC lidar has 
the advantages of long detection distance, high imaging 
speed and super-resolution imaging capabilities. 

The work was supported by the Hi- Tech Research and 
Development Program of China under Grant Project No. 
2011AA120101 and No. 2011AA120102. 
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